Do As antisites destroy the ferromagnetism of (Ga,Mn)As? 
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The effect of the inclusion of As antisites in the diluted magnetic semiconductor (Ga,Mn)As is 
studied within density functional theory in the local spin density approximation. In the case of 
homogeneous distribution of Mn ions we find that the ferromagnetism is largely weakened by the 
presence of the antisites. This is due to compensation of the free holes which mediate the long 
range ferromagnetic order. In contrast, when two Mn ions are coupled through only one As ion, 
ferromagnetic and antiferromagnetic order compete. In this case the magnetic ground state depends 
on: i) the position of the As antisites relative to the Mn, and ii) the As antisite concentration. We 
explain our results using a model of competing antiferromagnetic super-exchange and ferromagnetic 
double-exchange via localized Zener carriers. The strong dependence of the ferromagnetic order on 
the microscopic configuration accounts for the large variation in experimental data. 
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In the past ten years the study of diluted magnetic 
semiconductors (DMS) has been strongly revitalized 
since the discovery of ferromagnetic order in lui-j^Mn^^As 
and Gai_:rMn^As [|[ |, §. The potential utility of 
room temperature ferromagnetism in a semiconductor 
based system is enormous. On one hand, DMS could be 
replace the existing metallic magnetic elements in storage 
media Here their better compatibility with existing 
semiconductors and MBE-growth technology over their 
metallic counterparts is a great advantage. On the other 
hand DMS can be used as spin injectors into semiconduc- 
tors 1^ . This is a critical step in the physical realiza- 
tion of quantum computation based on the spin degree 
of freedom in a solid state device In fact, although it 
has been shown that spin can be coherently transported 
over several micrometers in GaAs ^, |l^, its injection 
from metallic magnetic contacts is largely unsuccessful if 
not impossible This is due to the large mismatch 

between the resistances of the metallic contacts and the 
semiconductor. The use of DMS solves this problem 0] 
and (Ga,Mn)As together with GaAs and (Al,Ga)As rep- 
resent to date the most promising material system for 
the injection, storage and manipulation of spins in semi- 
conductors. 

Although there is general agreement on the carrier- 
(hole-) mediated origin of the ferromagnetism in 
(Ga,Mn)As, the detailed mechanism is still debated. The 
Zener model in the mean field approximation [T^ pro- 
vides a good starting point. However its prediction of 
the Curie temperature as a function of Mn and hole 
concentration largely overestimates the actual value if 
a dynamic description of the Mn spins is considered 
[Q. Moreover recent density functional theory (DFT) 
results show that the p-d coupling in (Ga,Mn)As 
is very strong and the mean field approximation cannot 
be completely justified. Finally DFT calculations based 
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on the coherent-potential approximation (CPA) [jTsj sug- 
gest that the ferromagnetic order in (Ga,Mn)As may re- 
sult from the competition between ferromagnetic double- 
exchange and antiferromagnetic super-exchange interac- 
tions. 

Note that all these models agree on the following im- 
portant points: i) Mn^+ ions substitute the Ga'^+ cations 
in the zincblende lattice providing local 5 = 5/2 spins, 
and ii) there are free holes in the system although the 
actual concentration is much smaller than the density 
of Mn ions |j. This latter point suggests that some 
mechanism of compensation must be taking place, most 
likely the presence of intrinsic donor defects such as As 
antisites (Asca) or Ga vacancy-interstitial As pairs (Vca- 
Asi), which are usually found in low-temperature GaAs 
]l6[ . In the above models the compensation mechanism 
is implicitly incorporated in the mean field description of 
the GaAs valence band, since the hole concentration is a 
free parameter. As a result the local effects of magnetic 
and chemical disorder are neglected. 

In this paper we consider explicitly the effects of the in- 
clusion of AsGa in (Ga,Mn)As at different dilutions, and 
study how the chemical environment modifies the mag- 
netic interaction between the Mn ions. We perform DFT 
calculations within the local spin density approximation 
(LSDA), and consider large GaAs cells in which Mn ions 
and AsGa are inserted. We study the dependence of the 
magnetic coupling between Mn ions on: i) Mn concen- 
tration, ii) AsGa concentration, iii) relative positions of 
the Mn ions and the AsQa defects. The main result of 
our analysis is that, for a uniform distribution of Mn 
ions, AsGa antisites strongly weaken the ferromagnetic 
coupling. However, at least for moderate concentrations, 
the compensation does not follow the expected nominal 
valences of Mn and As (that is one AsGa does not com- 
pletely compensate the holes from two Mn ions). More- 
over when two Mn ions are separated by only one As 
atom (so that they occupy two corners of a zincblende 
tetrahedron with an As atom in the center), ferromag- 
netic coupling is possible even far above compensation 
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if As antisites are located at the other two tetrahedral 
positions to form Mn2As3 complex (see Fig.^). 

Our calculations are performed using the code siesta 
p7| , [T^ , , which is an efficient implementation of DFT- 
LSDA based on pseudopotentials and a numerical local- 
ized atomic orbital basis set. This method combines 
good accuracy and small computational cost compared 
to other methods based on plane- waves. We have de- 
scribed in a previous paper jlj] the calculation details 
and the prescriptions one needs to optimize the basis set 
and the pseudopotentials. Here we just mention that 
we use TrouUier-Martins pseudopotentials [ pOf with non 
linear core corrections pl[ and Kleinman-Bylander fac- 
torization and the Ceperley- Alder |23j form of the 
exchange-correlation potential. 

We construct 64 (cubic) and 32 (rectangular) atom 
GaAs cells in which we include two Mn ions and a vari- 
able number of Asca antisites. These correspond to Mn 
concentrations of x=0.0625 (to date the largest concen- 
tration obtained experimentally) and x=0.125 respec- 
tively. We investigate two different Mn configurations 
for each concentration: 1) the Mn ions occupy positions 
as far apart as possible (i.e. the corner and the middle 
of the cubic cell), 2) the Mn atoms occupy two corners 
of a tetrahedron and are coordinated through a single As 
ion (see Fig.||). We call these two configurations sepa- 
rated and close respectively. The cells are then periodi- 
cally repeated using 18 fc-points in the irreducible cubic 
Brillouin zone. Note that the periodic boundary con- 
ditions fix the magnetic coupling between Mn ions that 
occupy equivalent positions in adjacent cells to be ferro- 
magnetic but that the coupling between the two Mn ions 
in the same cell can be ferromagnetic or antiferromag- 
nctic. We calculate the ground-state properties of these 
systems for both ferromagnetic (FM) and antiferromag- 
netic (AF) alignment of the Mn ions within the cells. 

The energy differences between the AF and FM aligned 
configurations, ApA, and the magnetization per Mn ion, 
-^Mn, for the separated arrangement are presented in 
Fig.|l| as a function of the number of Asca antisites. The 
magnetization per Mn ion is defined to be half of the 
magnetization of the cell, calculated in the FM aligned 
phase. 

First we note that the ferromagnetic coupling is 
strongly weakened by Asca antisite doping. This is con- 
sistent with the picture of magnetic coupling mediated 
by free carriers (holes): AsQa antisites contribute elec- 
trons into the system and therefore compensate the holes. 
We also note that in the case of no antisites the fer- 
romagnetic order is much stronger in the case of large 
Mn concentration, which again supports the model of 
carrier-induced ferromagnetism. Within this model, the 
ferromagnetic coupling should disappear when the com- 
pensation is complete. However the figure suggests that 
the compensation mechanism does not follow the nominal 
atomic valence, since a single Asca antisite per cell is not 
sufficient to destroy the ferromagnetic coupling. Above 
compensation (one Asca for two Mn ions) antiferromag- 
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FIG. 1: (a) Energy difference between AF and FM align- 
ments, ApA, and (b) magnetization per Mn ion, Mmh, as a 
function of the number of Asca antisites in the cell: sepa- 
rated configuration. The horizontal line denotes the division 
between FM and AF alignment. 



netic coupling is obtained for large Mn concentration, 
while the system stays ferromagnetic at low concentra- 
tion. This is consistent with the onset of antiferromag- 
netic super-exchange coupling, the mechanism which is 
believed to be responsible for the magnetic order in the 
H-VI DMS ||2^, at compensation. Super-exchange is a 
short range interaction and therefore is less important in 
the low concentration limit where the Mn ions are well 
separated. It is worth noting that for a;=0.125 and large 
antisite concentrations {n AsGa=2) the magnetic order is 
quite sensitive to the actual position of the antisite re- 
spect to the Mn ions. Different AsQa arrangements may 
result either in ferromagnetic or antiferromagnetic cou- 
pling. 




FIG. 2: DOS for (Ga,Mn)As obtained for a 64 atom unit 
cell with two Mn ions and one Asca antisite: separated con- 
figuration and ferromagnetic alignment. The upper graph 
represents the majority and the lower the minority spins. 

If we now consider the magnetization per Mn ion 
(Fig.|l|b) we note that there is a monotonic increase of 
the total magnetic moment on antisite doping. This can 
be understood by looking at the density of states of Fig.|^ 
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(see also Ref.|^. At the Fermi energy the DOS of the 
majority spin band is the result of the hybridization be- 
tween the Mn d states with t2 symmetry and the p states 
coming mainly from Asca- In fact the two donor levels 
of isolated AsQa in GaAs are known to lie 0.54 eV and 
0.75 eV above the GaAs valence band [|l6) and to over- 
lap with the exchange split Mn d levels of low dilution 
(Ga,Mn)As 0. In contrast the DOS of the minority 
spin band is dominated by antisite levels, the first Mn d 
states lying far above Ep. Therefore by increasing the an- 
tisite concentration, Ep is shifted towards higher energy, 
thus increasing the occupation of the Mn d shell in the 
majority spin band. This enhances the total magnetic 
moment of the system. Finally note that a magnetiza- 
tion of 4.5/iB is in very good agreement with recent x-ray 
magnetic circular dichroism measurements psj . 

Let us now turn our attention to the close configu- 
ration. In this case we expect the short range super- 
exchange to play a more important role since the Mn 
ions are much closer to each other. Moreover the mag- 
netic coupling within a cell is expected to be quite sen- 
sitive to the position of the antisites with respect to 
the Mn ions. We consider the three different situa- 
tions sketched in Fig.^: a) the antisites are far from the 
Ga2Mn2Asi complex, b) one antisite occupies a tetrahe- 
dral site (GaiMn2As2), c) two antisites occupy the tetra- 
hedral sites (Mn2As3). In Fig.^ we show ApA and Mun 
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FIG. 3: Close Mn configuration. The three antisite arrange- 
ments considered in the text: a) Ga2Mn2Asi, b) GaiMn2As2, 
c) Mn2As3. 




(b) 




0.5 1 1.5 2 2.5 3 

n As^ 



FIG. 4: (a) Energy difi'erence between the AF and FM align- 
ments Afa, and (b) magnetization per Mn ion Mmu as a 
function of the number of Asca antisites in the cell: Close 
Mn arrangement. The symbols •, * and o represent arrange- 
ments (a), (b) and (c) of Fig.^ respectively. 



we propose that super- exchange coupling stabilizes the 
AF phase as soon as the free-holes are completely com- 
pensated. However this mechanism is extremely short 
range and therefore effective only when the Mn ions are 
separated by one As. 

We now turn our attention to situations (b) and (c). In 
both cases the ferromagnetic alignment is stable and al- 
most insensitive to the total Asoa concentration. This is 
a strong indicator that the dominant interaction in these 
cases is one that depends only on the local chemical prop- 
erties. To shed some light on this aspect we present the 
results of Miilliken population analyses of the 

system. The Miilliken orbital (or atomic) population is 
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as a function of the AsQa concentration for the different 
antisite arrangements of Fig.^. For the sake of brevity 
we present data only for a;=0.0625, noting that similar 
conclusions can be made also for a;=0.125. From FigJJ it 
is very clear that, when the Asca defects are in the vicin- 
ity of the Mn ions (configurations (b) and (c) of Fig.||), 
the system is almost insensitive to the Asca concentra- 
tion. In contrast when the Asca antisites are far from 
the Mn ions the magnetic coupling undergoes a FM to 
AF transition with increasing antisite concentration. 

Consider first configuration (a). In this case both the 
ApA and A/mii curves look very similar to the curves we 
found for the separated arrangement (Fig.^. However, 
in stark contrast with the separated case, an AF align- 
ment can now also be found for large Asca concentra- 
tions. Recalling the fact that the Mn d shell is more 
than half-filled, and that it is antiferromagnetically cou- 
pled with the As p shell of the intermediate atom p| . 



TABLE I: Miilliken atomic and orbital populations for the 
Mn ions and the intermediate As atom of the complexes of 
Fig.^. The Mn concentration is a;=0.0625. The symbols t 
and I correspond to majority and minority spin respectively. 
The populations are in units of the electronic charge |e|. 

the projection of the charge density onto a particular or- 
bital (or atom). Although the absolute values depend 
on the chosen basis set, it is a useful quantity for under- 
standingthe charge distribution within a given material. 
In table || we present the Miilliken orbital population for 
the two Mn ions and the intermediate As ion of the com- 
plexes of Fig.^ The most obvious feature is the larger 
polarization of the As-p orbitals in Mn2As3 compared 
with Ga2Mn2Asi. This is accompanied by a small de- 
crease in Mn-d polarization, although this latter effect 
could be an artifact from the overlap component of the 
orbital population |l3] . It is also interesting to note that 
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the total atomic population of the As atom increases go- 
ing from (a) to (b) to (c), while the p component of the 
population decreases. Therefore if we start from the sit- 
uation in which two Asca antisites are located far from 
the Mn-As-Mn complex, then we move each antisite in 
turn to one of the two other corners of the tetrahedron, 
then i) the charge on the middle As atom increases, ii) 
the spin-polarization of the p shell of the middle As atom 
increases, and iii) the total population of the p-shell of 
the middle As atom decreases. And most importantly 
the magnetic coupling changes from antiferromagnetic to 
ferromagnetic. 

In 1960 de Gennes observed |^ that in an antiferro- 
magnetic crystal the presence of a bound carrier (elec- 
tron or hole) which is Zener coupled to the local spins 
always induces a distortion in the antiferromagnetic lat- 
tice. In that spirit we propose that the observed tran- 
sition from antiferromagnetic to ferromagnetic coupling 
between Ga2Mn2Asi and Mn2As3 results from the on- 
set of ferromagnetic double-exchange coupling mediated 
by a bound Zener carrier. This is further supported by 
recent calculations in which we show that the ac- 
tual ground state for GaiMn2As2 is not one of perfect 
ferromagnetic alignment, but instead the Mn magnetic 



moments are canted with respect to each other. 

In conclusion we have shown that the Asca antisites' 
inclusion in (Ga,Mn)As can result in a variety of different 
behaviors depending on the microscopic arrangement of 
the Mn ions and the Asca antisites. In particular we have 
shown that when the Mn ions are uniformly distributed 
in the crystal, Asca antisites weaken the ferromagnetic 
order. In contrast if the Mn ions occupy two corners 
of the zincblende tetrahedron several magnetic arrange- 
ments are possible depending on the positions of the an- 
tisites. In particular ferromagnetic coupling is obtained 
if antisites occupy the other positions in the tetrahedron. 
This suggests that a way to obtain high Tc in (Ga,Mn)As 
is either to reduce the As antisite concentration, or to de- 
posit Mn and As antisites with a strongly inhomogeneous 
distribution. 
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